FACILITY FORM 602 


NASA TECHNICAL NOTE 


Os 

Os 

sO 

CO 



J 



GPO PRICE $ 

CFSTI PRICE(S) $ 


#62s, ON UJ 8 8 


ST’ _ 

rm^ui 

(PAGES) 

(CODE) 

(NASA CR OR TMX OR AD NUMBER) 

fCATEGORY) 


Hard copy (HC) _ 

Microfiche (MF) & 

ff 653 July 86 


ANALYSIS OF USE OF JUPITER GRAVITY 
TURNS FOR TAILORING TRAJECTORIES 
TO SATURN WITH CONSIDERATION 
OF LAUNCH OPPORTUNITIES 


by Edwin F. Harrison and Charles H. McLellan 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • NOVEMBER 1966 


NASA TN D-3699 


ANALYSIS OF USE OF JUPITER GRAVITY TURNS FOR TAILORING 

TRAJECTORIES TO SATURN WITH CONSIDERATION 

OF LAUNCH OPPORTUNITIES 

By Edwin F. Harrison and Charles H. McLellan 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.50 


ANALYSIS OF USE OF JUPITER GRAVITY TURNS FOR TAILORING 


TRAJECTORIES TO SATURN WITH CONSIDERATION 
OF LAUNCH OPPORTUNITIES 

By Edwin F. Harrison and Charles H. McLellan 
Langley Research Center 

SUMMARY 

A trajectory study has been conducted to determine several launch opportunities for 
Earth to Saturn missions, especially those times when the intermediate planet Jupiter can 
be passed closely enough to use its strong gravitational field to make favorable, significant 
changes in the interplanetary-trajectory conditions. The effect of hyperbolic-excess 
velocity and flight -path angle at departure from Earth on both launch date and trip time 
are studied for the gravity-turn trajectories. These trajectory characteristics are com- 
pared with direct Earth to Saturn trajectory requirements. Passage distance, turning 
angle, and velocity gained at Jupiter for the gravity-assisted trajectories are examined. 
This study is based mainly on coplanar trajectories; however, an evaluation of noncoplanar 
effects is included. 

Jupiter gravity turns can be utilized for missions to Saturn with launch dates in 
1976, 1977, 1978, and 1979, and again starting in about 1996. With the approximate 
minimum-energy requirement for direct missions to Saturn, a reduction in trip time of 
about 40 percent was realized for the gravity-turn trajectories over the direct trajec- 
tories. By employing a gravity-turn mode, missions to Saturn can also be made with 
velocity levels at about 1 km/sec lower than the minimum value required for the direct 
trajectories. This reduction in required velocity can double the payload fraction for a 
fly-by Saturn mission if a launch vehicle with approximately the same capability as the 
Saturn V is used. 


INTRODUCTION 

As planetary exploration advances to more distant planets, the trip times become 
increasingly long and the energy requirements large. For example, a minimum -energy 
trajectory from Earth to Saturn requires about 6 years and a velocity of about 10.4 km/sec 
at departure from the Earth's sphere of influence. Shorter trip times can be obtained at 
the expense of an increased energy requirement, which, of course, is very costly in terms 
of payload fraction of the launch vehicle. Another approach that looks attractive for some 



of these missions (refs. 1 and 2) is to utilize the free energy that can be derived from the 
gravitational field of the intermediate planet Jupiter when it is possible to pass near this * 
planet before reaching Saturn. This gravitational effect is referred to in the literature 
as gravity turn, gravity assist, planet swingby, or planetary attraction. These studies of 
Earth-Jupiter-Saturn trajectories, as well as studies of Earth-Venus-Mars combinations 
(refs. 3, 4, and 5), have shown that the use of gravity-assisted trajectories can offer 
advantages, such as reductions in required launch velocity, trip time, or arrival velocity, 
when compared with the values for direct missions. If a huge planet like Jupiter can be 
used, its strong gravitational field can be very effective in tailoring interplanetary 
trajectories. 

The use of a gravity turn at Jupiter for an Earth-Saturn mission requires particular 
combinations of the relative positions of all three planets; therefore, the number and fre- 
quency of the launch opportunities are limited. There has been very little consideration 
given to launch date in the literature * to disclose when full benefits from Jupiter gravity 
turns can be realized. However, one launch date for a specific launch velocity is dis- 
cussed in reference 2. 

The purpose of the present study is to determine, for Jupiter gravity-assisted tra- 
jectories to Saturn, launch dates for several launch opportunities and the corresponding 
trip times for a range of hyperbolic-excess velocities and to compare the results with 
direct-mission requirements. The passage distance, turning angle, and velocity gained 
at Jupiter for the gravity-assisted trajectories are also presented. 

ANALYSIS 

Coplanar Trajectories 

A two-dimensional, two-body analysis was performed for the main part of the pres- 
ent study in order to determine approximate launch opportunities and the associated 
characteristics of Jupiter gravity-assisted trajectories to Saturn. In order to analyze the 
trajectories, they were divided into two main parts: (1) the departure from Earth to 

arrival at Jupiter and (2) the gravity turn and departure from Jupiter to arrival at Saturn. 
The geometry involved in these heliocentric trajectories is presented in figure 1. The 
planets are assumed to be point masses with their spheres of influence shrunk to zero for 
the heliocentric trajectories, and thus the trajectories go from the center of one planet to 
the center of another. However, for the gravity-turn maneuver at Jupiter, consideration 
must be given to the gravitational field of the planet and passage distance, as shown in 
figure 2. 

1 After submission of this report for publication, the authors learned of a similar, 
concurrent study by Jerry M. Deerwester entitled "Jupiter Swingby Missions to the 
Outer Planets," AIAA Paper No. 66-536, June 1966. 
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In general, the approach followed in this study is to simplify the analysis without 
seriously degrading its accuracy. Since the Earth-Jupiter central angles to be considered 
in this analysis will not be near 180°, the out-of-the-ecliptic inclination angles and veloc- 
ity requirements are small and can be neglected for first-order approximations with good 
accuracy for mission feasibility studies; hence a coplanar trajectory can be used. 

Since the eccentricity of Earth is very small, its orbit can be assumed to be cir- 
cular with little error being introduced into the interplanetary trajectories. However, 
the eccentricities of Jupiter and Saturn are appreciable and must be accounted for. 

The symbols used in the trajectory equations are defined in appendix A, and the 
equations used for the coplanar trajectories are given in appendix B. In order to evaluate 
the effect on the analysis of assuming that the trajectories are coplanar, some calculations 
were made for noncoplanar trajectories; the equations are presented in appendix C. 

Earth to Jupiter trajectories .- The coplanar trajectory considered for this study 
begins at Earth, where the hyperbolic -excess velocity of the spacecraft and its helio- 
centric flight-path angle are specified. This velocity and angle are then used in the equa- 
tions given in appendix B to calculate trajectories to Jupiter's elliptical orbit, where the 
arrival heliocentric radius depends on launch date and trip time. In order to determine 
the matching conditions between these trajectory parameters, a graphic iteration was 
employed. In this procedure, calculations of trip time, central angle, and heliocentric 
radius from the trajectory are uniquely matched with launch date and other comparable 
ephemeris data (refs. 6 and 7) of Earth and Jupiter. A detailed explanation of this 
matching procedure along with a graphic illustration is given in appendix B. 

Jupiter to Saturn trajectories .- As the spacecraft passes Jupiter, the gravitational 
field of the planet causes changes in heliocentric direction and velocity, as schematically 
illustrated in figure 1. The change in the heliocentric velocity occurs because the arrival 
and departure hyperbolic-excess velocity vectors are different in direction but equal 
in magnitude. The spacecraft passes on the back side of the planet in an effort to increase 
the heliocentric velocity. The closest passage distance (see fig. 2) has to be determined 
such that the spacecraft is located on an encounter trajectory with Saturn, the target 
planet. 

Because of the number of variables involved, both calculations and graphs are used 
to determine the proper passage distance at Jupiter. Trajectories are computed from 
the equations in appendix B by using the determined arrival conditions and various selected 
values for the passage distance and the heliocentric radius to Saturn. These trajectory 
calculations are then plotted along with Jupiter and Saturn ephemeris data, and matched 
conditions are obtained. (See appendix B.) The matched trajectory conditions include: 
trip time, central angle, heliocentric radius, and passage distance. With these quantities 


3 



determined, the desired coplanar trajectory to Saturn involving a gravity assist at Jupiter 
is established. 

It should be noted that this two-dimensional analysis provides only approximate 
conditions for the gravity-assisted trajectories. For the purpose of determining the accu- 
racy of this analysis a noncoplanar or a simplified three-dimensional analysis was also 
performed. 


Noncoplanar Trajectories 

When the orbital inclination of Jupiter is taken into account, the orbital plane of the 
spacecraft, which is assumed to lie in the ecliptic for the two-dimensional analysis, has 
to be changed in order to intercept Jupiter. This orbital plane change requires additional 
velocity, which is considered for convenience to be applied either at Earth upon departure 
or at Jupiter's ascending node. At either location, all the coplanar characteristics 
between Earth and Jupiter, except of course the total velocity requirement, can be main- 
tained by rotating the velocity vector of the spacecraft coplanar trajectory in a plane 
normal to the elliptic. These plane changes are shown graphically in figure 3 and the 
equations are presented in appendix C. 

If the orbital plane of the spacecraft is inclined to the ecliptic upon departure from 
Earth, then it will intercept Jupiter's orbital plane upon arrival at the angle of declination 
\2 j. (See fig. 4.) Hence, the equations for the coplanar trajectory become inapplicable 
at Jupiter and for the Jupiter to Saturn phase of the mission. Equations were developed 
(see appendix C) to allow for the angles between the approach and departure planes of the 
spacecraft at Jupiter. Noncoplanar trajectories were then calculated from these equa- 
tions and matched with planetary ephemeris data, in a manner similar to that used for the 
coplanar trajectory between Jupiter and Saturn. 

RESULTS AND DISCUSSION 
General Effects of Gravity Turns 

Before considering specific launch opportunities and trajectories, some of the gen- 
eral characteristics of the Jupiter gravity-turn mode for circular planetary orbits are 
examined. The closest passage distance of the spacecraft at Jupiter is one of the prime 
factors because it has a strong effect on the Jovian departure flight- path angle, which 
results in heliocentric velocity gains, as shown in figures 5 and 6. For each value of the 
Earth departure hyperbolic-excess velocity Vh j considered, the maximum heliocentric 
velocity gained at Jupiter occurs when the spacecraft passes near the surface of the 
planet. This spacecraft position corresponds to a Jovian departure flight-path angle of 0°. 
It might be first thought that these conditions would yield the minimum trip time to Saturn. 
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However, the data in figure 7 show that the minimum trip times usually occur when the 
passage distance is several radii from the planet, where the velocity gained is less than 
the maximum possible value. Thus, as the passage distance increases to some extent, 
the associated change in Jovian departure flight-path angle can more than compensate for 
the decrease in velocity gained. 

A comparison of the minimum trip times in figure 7 with direct trip times is pre- 
sented in figure 8. A large reduction in trip time for a specified value of V^i or a 

decrease in the required value of for a selected trip time is shown for the gravity- 

turn trajectory as compared with the direct trajectory. In order to use such gravity 
turns, however, the ephemeris of the planets must be considered and launch opportunities 
determined. 


Position Requirements of Planets 

The relative positions of the planets for a 30-year period and the approximate angle 
requirements of them at launch time for the Jupiter gravity-turn trajectories are shown 
in figure 9. The ephemeris angles presented are the differences in the heliocentric longi- 
tudes of Jupiter and Earth along with the corresponding differences in the positions of 
Jupiter and Saturn. (See refs. 6 and 7.) The Earth-Jupiter synodical cycle is about 
1 year, whereas the Jupiter -Saturn synodical period is nearly 20 years. The relative 
locations of the planets required for the gravity-turn trajectories were determined as 
follows: 

Trajectory transfer angles were calculated for Earth departure hyperbolic-excess 
velocities ranging from 10.363 km/sec (34 000 fps), representing the approximate 
minimum-velocity requirements for a direct trajectory to Saturn, to an arbitrary value of 
12.192 km/sec (40 000 fps); the departure heliocentric flight -path angle at Earth was 
initially 0°. These computed angles were then transferred, by using the motion of planets, 
to the angles required between the planets at launch. The shaded areas in figure 9 repre- 
sent the required angles. When the ephemeris data lines for the planets pass through 
the required-angle regions for the trajectories, launch opportunities become available. 

The Jupiter-Sun-Earth angle requirement can be met almost every year; however, there 
is only about a 3-year period every 20 years in which the required and available Saturn- 
Sun- Jupiter angles can be matched. From this figure it can be seen, therefore, that launch 
opportunities for Jupiter gravity-assisted trajectories will be available in 1976, 1977, 

1978, and 1979, and again starting in 1996 for the velocities considered. 
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Effects of Hyperbolic-Excess Velocity on Trajectory Characteristics » 

For the 1978 launch opportunity, the effect of the hyperbolic-excess velocity on 
other trajectory characteristics such as launch window, turning angle and closest passage 
distance at Jupiter, heliocentric profiles, and trip times for Jupiter gravity trajectories 
to Saturn are presented in figures 10 to 12. These results are based on the coplanar tra- 
jectory analysis. 

The data of figure 10 show that the range of hyperbolic -excess velocities considered 
for tangential launches from Earth’s orbit provides a launch window of about 8 days in 
1978 for both Jupiter gravity-turn and comparable direct trajectories to Saturn. Although 
only one launch date is shown to be available for a given value of V^l because yj = 0°, 

the launch window can be enlarged by varying y^. For the present launch conditions, the 
corresponding Earth-Saturn trip times are shown in the upper part of this figure. A vari- 
ation of the hyperbolic -excess velocity at Earth was shown to have the greatest effect on 
trip time at the lower velocities. At the approximate minimum velocity for a direct 
mission, the trip time obtained for the gravity-turn trajectory is as much as 40 percent 
lower than that for the direct trajectory. When both types of trajectories have the same 
trip time, say 5 years, an important reduction of about 1 km/sec or approximately 10 per- 
cent in hyperbolic-excess velocity is indicated for the gravity-turn trajectory. Also pre- 
sented in figure 10 is the injection velocity, which is equal to the square root of the sum 
of squares of the hyperbolic-excess velocity and the escape velocity. Thus, the 10- 
percent reduction in hyperbolic-excess velocity corresponds to lowering the injection 
velocity by about 5 percent. 

Another possible use of the gravity turn is to provide additional launch windows, 
which would be approximately a month ahead of the direct-mission launches. For 
example, if a value of 11 km/sec can be attained for V h j, it would be possible, in 1978, 
to launch on October 10 by using a gravity-turn trajectory or on November 5 by using a 
direct trajectory to Saturn. (A shift of a few days in these launch dates will, of course, 
occur when an exact analysis is used.) 

The closest passage distance and turning angle which correspond to the trajectory 
characteristics of figure 10 are presented in figure 11. The hyperbolic -excess velocity 
is shown to have a pronounced effect on the passage distance required at Jupiter. This 
effect was strongest at the lower velocities where a decrease of 1 km/sec caused an 
increase in the passage distance of about 30 Jupiter radii, while the corresponding turning 
angle decreased only about 8°. 

Illustrated in figure 12 is the effect of the hyperbolic -excess velocity on the helio- 
centric profile of the gravity-turn trajectory from Earth to Saturn. No appreciable effects 
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on the appearance of the trajectories are evident, even though there is nearly a year 
difference in the trip times. 


* 


Comparison of Coplanar and Noncoplanar Trajectories 

All the preceding results are based on coplanar trajectories determined from a two- 
dimensional analysis. In table I, a comparison is made of some of the results obtained 
from this analysis with those determined by a similar analysis in reference 2 for a partic- 
ular launch. The present data agree well with reference results. The small difference 
in them can generally be contributed to the difference in heliocentric radii used. 

In order to determine the uncertainties involved in the results of the coplanar tra- 
jectories, they are compared with noncoplanar results in table II. The coplanar results 
are in good agreement with those of the noncoplanar trajectories; therefore, the coplanar 
results presented in this report can be relied on to represent good approximations and 
proper trends for mission planning. A comparison is not made of the present results 
with the Earth-Jupiter-Saturn trajectory results given in reference 8, because a 30° error 
was made in determining the heliocentric position of Saturn in that report, as noted in 
reference 9. 

Whereas the departure velocity from the Earth's orbit is assumed to be the same 
for the coplanar and noncoplanar trajectories, a velocity increment is required to change 
the trajectory plane so that the trajectory will intercept Jupiter. As discussed in appen- 
dix C, this velocity increment can be added in such a way that the plane change is obtained 
without a change in the heliocentric velocity or the flight-path angle of the spacecraft. 

For convenience, the plane change is considered to occur either at the departure from 
Earth or at Jupiter's node. In figure 13, the additional velocity required for these plane 
changes is presented as a function of launch date for a given hyperbolic-excess velocity. 
For the gravity-turn trajectories, the additional velocity required is very small when it is 
applied at departure from Earth, because here it can be combined with the coplanar 
hyperbolic-excess velocity and used as a single impulsive velocity. At Jupiter's node the 
velocity increment shown represents the total propulsion required there. Regardless of 
the plane-change location, the comparable direct trajectories shown require a larger 
velocity increment than the gravity-turn trajectories, mainly because they have transfer 
angles of almost 180°. However, none of the velocity increments given in figure 13 are 
included in the other results presented in this report, because they can be included after 
the trajectories are determined.. 

Both the angular approach and the departure inclination angle at Jupiter affect the 
Jupiter-Saturn leg of the gravity-turn trajectory. The approach characteristics depend 
on the location of the spacecraft when the plane change is initiated before reaching Jupiter. 
Because of the small velocity increment required for a plane change at departure from 
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Earth for the gravity-turn trajectories (see fig. 13), this location was used for the non- 
coplanar trajectories. The required values of the parameters at departure, however, can * 
be attained without propulsive velocity by making proper use of the gravity-turn maneuver , 
even when Saturn's orbital inclination is included. For example, an upward deflection of 
the departure velocity can be obtained by approaching Jupiter below its orbital plane. 

Effects of Launch Year on Trajectory Characteristics 

As indicated in figure 9, launch opportunities for gravity-assisted trajectories to 
Saturn will be available in 1976, 1977, 1978, 1979, and again in about 20 years. The effect 
of these opportunities on trip time is presented in figure 14 along with comparable direct 
trip times required for the same hyperbolic-excess velocity. For these launches with 
j of 10.363 km/sec, the trip time for any particular gravity -turn trajectory is some 
35 to 45 percent less than that of the direct trajectory. The shortest trip times for the 
gravity-turn trajectories are in 1977 and 1978. Also, there is good agreement between 
the trip times of the coplanar and the noncoplanar trajectories. 

The corresponding passage distance and turning angle at Jupiter required for the 
gravity-assisted trajectories shown in figure 14 are a strong function of launch year, as 
illustrated in figure 15. The closest passage distance occurs for the 1976 launch date; 
for each succeeding launch date, the passage distance increases very rapidly. No launch 
opportunities in which the gravity turn could be used were found prior to 1976, even 
though a range of values of \ were considered. For the 1979 launch date shown in 
figure 15, the passage distance is over 100 Jupiter radii and the turning angle is only about 
10°. These results indicate that for launches later than 1979 there probably would be 
only small advantages gained (i.e., reduced trip time) by using a gravity turn, unless the 
hyperbolic-excess velocity at Earth is increased. (See fig. 11.) The turning angle and 
passage distance for the coplanar trajectories shown in figure 15 agree reasonably well 
with those of the noncoplanar trajectories, except for the passage distance at the 1979 
launch date. For this date, the passage distance is so far from the planet that its effect 
on the trajectory would be rather small anyway. 

Presented in figure 16 are the actual velocity gained and the angle change at Jupiter 
associated with the launches from 1976 to 1979. These trajectory changes are strongly 
influenced by the passage distance, and their maximum values are obtained for the 1976 
launch date, which had the closest passage distance to Jupiter (fig. 15). For this launch 
date the magnitudes of these parameters are large: about 11.5 km/sec is gained in space- 
craft velocity and a change of approximately 63° occurs in the flight -path angle. For this 
trajectory almost maximum utilization is made of this free velocity which is available at 
Jupiter. (See fig. 5.) 
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As the velocity gained at Jupiter increases, the arrival hyperbolic-excess velocity 
* at Saturn becomes larger as shown in figure 17. A comparison of these arrival velocities 
with those of direct trajectories for comparable trip times has not been made here. How- 
ever, such a comparison has been made in reference 10 and, although the details are not 
shown, the comparison indicates that arrival velocities are so high that if an orbiter mis- 
sion to Saturn is desired, the Jupiter gravity turns become less attractive than the direct 
missions. 

In figure 18, heliocentric profiles of the spacecraft trajectory for the launches from 
1976 to 1979 for Earth -Jupiter -Saturn missions are presented. The Earth-Jupiter central 
angles of the trajectories have essentially the same value for all the launches, even though 
they start at different points. However, the Jupiter-Saturn central angle decreases con- 
siderably for each succeeding launch year because the trip time remains about the same 
(fig. 14), while the gain in velocity and the flight -path angle change at Jupiter (fig. 16) 
decrease. 


Solar Occultation 

As the spacecraft approaches Saturn, the line of sight to Earth can possibly be 
occulted by the Sun. The effect of trip time on this solar occultation is presented in fig- 
ure 19. The trip time should be something other than about 3.1 years for the launch con- 
ditions considered in order to avoid solar occultation. In order to insure a clear line of 
sight, the Saturn-Sun-Earth angle should probably be less than about 170° because of the 
extension of the Sun's corona. 


Variation of Flight-Path Angle 

All the preceding results have been for a tangential launch (i.e., yj = 0°), and thus 
the launch window is very narrow for a particular year. In order to widen this launch 
window, the departure flight -path angle at Earth was varied for a constant hyperbolic- 
excess velocity. The results for the gravity-assisted trajectories are shown in figure 20. 
A variation from -4° to 4° in the flight-path angle produced a launch window of about 15 to 
20 days, depending on the launch year. This extension of the launch window should provide 
enough time for launch considerations. 

The effect of flight-path angle at Earth on trip time is presented in figure 21. Even 
though the hyperbolic-excess velocity was held constant, the departure heliocentric veloc- 
ity at Earth of the spacecraft decreased as the flight-path angle changed. Therefore, the 
trip time increased as the flight -path angle varied from 0°. As pointed out earlier, the 
minimum trip times are for the launches in 1977 and 1978. 

The combination effect of hyperbolic-excess velocity and flight -path angle on the 
launch window is illustrated in figure 22 for the launch year 1978. As the flight -path angle 
varies from -4° to 4°, the number of available days for launch is shown to be about 19 and 
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is almost independent of the hyperbolic-excess velocity considered. For a particular 
flight -path angle, the launch date is moved forward only about 3.5 days for a 1 km/sec ♦ 
increase in velocity. 

Figure 23 shows the effect of hyperbolic-excess velocity and flight -path angle on 
trip time. As expected, the shortest trip times occur when the hyperbolic-excess velocity 
increases and the flight -path angle approaches 0°. 

Payload Comparison 

A comparison of the payload fraction for a direct mission with that for a Jupiter 
gravity-turn mission to Saturn is presented in figure 24. The line shown represents the 
capability of a launch vehicle which is approximately the same as that of a Saturn V. The 
circular- and square-shape symbols on this line represent the minimum velocity require- 
ments for a Jupiter gravity-turn and a direct trajectory to Saturn, respectively, for a 
typical launch opportunity. These results indicate that a payload fraction for a fly-by 
Saturn mission can be doubled by employing a Jupiter gravity-turn trajectory. Further- 
more, the trip time would be shorter for the gravity-turn trajectory, according to fig- 
ure 10. 


Other Considerations 

This study has been concerned with only pure gravity turns at Jupiter, which have 
been demonstrated to be very advantageous for Earth to Saturn trajectories. However, 
more improvement in the performance characteristics of these trajectories might be 
obtained in future studies if the gravity turn is augmented with propulsion. In refer- 
ence 11, Mars nonstop round-trip trajectories with gravity and propulsion gravity turns 
are analyzed. Compared with gravity turns, propulsion-augmented gravity turns near 
Mars markedly expanded the range of possible trip times and favorable launch dates and, 
in some cases, reduced the required total propulsion velocity increment. 

Another aspect of the Jupiter gravity-turn trajectory that appears to need future 
consideration is the guidance and communication requirement for these missions. One 
study (ref. 12) has been made for an Earth-Venus-Mercury mission, and no areas of 
serious difficulties with regard to tracking, communication, or guidance were found. The 
results, however, indicated that combinations of some of the other planets would be more 
sensitive to errors in guidance. 


CONCLUDING REMARKS 

A two-dimensional analysis has been conducted on the utilization of the gravitational 
field of Jupiter for tailoring trajectories to Saturn. For these missions, launch dates in 
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, 1976, 1977, 1978, and 1979 have been determined for a value of the departure hyperbolic- 
excess velocity at Earth of near the minimum value required for a direct mission to 
Saturn. A reduction in the Earth-Saturn trip time of approximately 40 percent was 
realized for these gravity-assisted trajectories over the direct trajectories. The closest 
passage distance to Jupiter varied from very near the surface of the planet to more than 
100 Jupiter radii as the launch opportunity advanced from 1976 to 1979. 

For the gravity-turn analyses, the results obtained from the coplanar -trajectory 
equations agreed well with the noncoplanar values. 

A variation of the hyperbolic-excess velocity at Earth was shown to have the great- 
est effect on trip time at the lower velocities. Gravity-assisted trajectories to Saturn 
can actually be made with values of the hyperbolic-excess velocity which are a significant 
10 percent lower than the minimum velocity required for a direct trajectory. This reduc- 
tion in required velocity can double the payload fraction for a fly-by Saturn mission if a 
launch vehicle having approximately the same capability as a Saturn V is used. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 22, 1966. 
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APPENDIX A 


SYMBOLS 

a semimajor axis, kilometers 

E eccentric anomaly, radians 

e eccentricity 

F hyperbolic anomaly, radians 

I g p specific impulse, meters/second 

ijE orbital inclination of Jupiter to ecliptic, degrees 

ij g inclination angle between orbital planes of spacecraft and Earth at departure, 

degrees 

ig j declination angle between orbital planes of spacecraft and Jupiter at Jovian 

arrival, degrees 

igj inclination angle between orbital planes of spacecraft and Jupiter at Jovian 

departure, degrees 

l heliocentric longitude of planet, degrees 

1 E ,1 heliocentric longitude of Earth at time of departure from it, degrees 

l j 3 heliocentric longitude of Jupiter at time of departure from it, degrees 

Q velocity vector connecting arrival with departure hyperbolic-excess velocity 

vectors at Jupiter (see fig. 4), kilometers/second 

Rj radius of Jupiter, kilometers 

Rp periapsis distance of spacecraft from center of Jupiter, kilometers 

r radial distance of body from center of Sun, kilometers 
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APPENDIX A 

r a apoapsis distance, kilometers 

rp periapsis distance, kilometers 

T date with respect to orbital position of planet, Julian date 

T e i date at departure from Earth, Julian date 

Tj 3 date at departure from Jupiter, Julian date 

t time from periapsis, years 

t^ 2 trip time from Earth to Jupiter, years 

tg 4 trip time from Jupiter to Saturn, years 

V velocity of body relative to Sun, kilometer s/second 

V c circular velocity of Jupiter at surface of planet, kilometer s/second 

V h hyperbolic-excess velocity relative to planet, kilometer s/second 

AVh i velocity increment required to be combined with coplanar hyperbolic excess 
at Earth to make a plane change, kilometer s/second 

V n velocity of spacecraft at Jupiter's node, kilometers/second 

AV n velocity increment required at Jupiter's node to make a plane change, 

kilometers/second 

X axis which lies along intersection of orbital planes of spacecraft and Jupiter 

(see fig. 4 ) 

x direction of vernal equinox 

Y axis which is perpendicular to X-axis and lies in Jupiter's orbital plane 

Z axis which is perpendicular to Jupiter's orbital plane 
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APPENDIX A 


z distance out of ecliptic, kilometers 

z' distance between planes of Jupiter and Saturn upon arrival at Saturn, 

kilometers 

a gravity -turn angle with respect to Jupiter, degrees 

/3 angle between hyperbolic -excess velocity vector of spacecraft and orbital 

velocity vector of planet, degrees 

y flight -path angle relative to local normal of heliocentric radius vector, 

degrees 

y n flight -path angle of spacecraft at Jupiter's node, degrees 

9 orbital central angle measured from periapsis, degrees 

2 Earth-Jupiter central angle, degrees 

0g 4 Jupiter -Saturn central angle, degrees 

9-, central angle between Earth at departure and ascending node of Jupiter, 

i ,n 

degrees 

/i gravitational constant of Sun, kilometers ^/second^ 

celestial latitude of Jupiter at arrival, degrees 

a ratio of structural weight to booster stage weight 

angle between hyperbolic-excess velocity vector of spacecraft and Jupiter's 
plane, degrees 

SI longitude of ascending node, degrees 

Subscripts: 

1 trajectory conditions of spacecraft at departure from Earth 

2 trajectory conditions of spacecraft at arrival at Jupiter 
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3 trajectory conditions of spacecraft at departure from Jupiter 

4 trajectory conditions of spacecraft at arrival at Saturn 

E Earth's orbital conditions 

J Jupiter's orbital conditions 

5 Saturn's orbital conditions 
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APPENDIX B 


EQUATIONS USED FOR THE COPLAN AR TRAJECTORIES 

The type of coplanar trajectory considered in this paper may be divided into two 
parts: (1) the departure from Earth to arrival at Jupiter and (2) the gravity turn and 
departure from Jupiter to arrival at Saturn. Most of the equations required to define the 
characteristics of this trajectory have been developed (ref. 13) or depend only on simple 
geometric laws; hence, they will be presented herein with a minimum of description. 

Earth to Jupiter Trajectories 

The spacecraft is assumed to depart from Earth on an interplanetary trajectory to 
Jupiter (fig. 1) with a specified hyperbolic-excess velocity V h j and flight-path angle 
yj. The radial distance from the Sun at the departure point r j is assumed to be 
1.49 x 10® km, but the elliptical heliocentric radius of Jupiter *2 is a variable. By 
employing all these conditions, trajectories can be computed from the following equations. 
These equations are arranged in the order required for a straightforward, step-by-step 
calculation procedure. 

The departure heliocentric velocity at Earth, based on the geometry in figure 1 and 
the laws of cosines and sines, is 

V - V + v 2 - 2V E V M cos h < B1 > 

where 


/3j = 180 - yj - sin 


-1 


V V M 


sm 


y i 


The next seven equations are some of the fundamental elliptic orbit element relations 
as given in reference 13. 

Semimajor axis: 


a = 


r l 


2 - 


r l V l 2 

M 


(B2) 


Eccentricity: 


= 1 - ^-(2a - rj) cos^yj 
a 2 


(B3) 
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Perisapsis distance: 


Apoapsis distance: 


Time from periapsis: 


r p = a(l - e) 


(B4) 


r a = a(l + e) 


(B5) 


Orbital central angle: 


Eccentricity anomaly: 


6 = cos' 


•1 - ae^ - rj 


E = cos 




(B6) 


(B7) 


t = (E - e sin E) 


(B8) 


The following two equations are employed to transfer the reference point from the 
periapsis for the central angle and travel time to the departure point at Earth: 


Earth-Jupiter central angle: 


e l,2‘ e 2* e l ( ±6 1 ^ ^ (B9) 

Earth-Jupiter trip time: 

ti >2 = t 2 ±t l (^1 f ° r Ty l) ( Bl °) 

It should be noted that and tj are determined by using r = r^ = 1.49 x 10® km and 
6 2 and t 2 are determined by using r = in equations (B6) and (B8). 

In order to determine the launch date and the corresponding trip time to Jupiter for 
the particular trajectory considered (i.e., the values of V h j and y^ are specified) the 
calculated results from equations (B9) and (BIO) are plotted in sketch (a) along with some 
planetary ephemeris data taken from references 6 and 7. The resultant angle l j - Z E i 
is the difference between the heliocentric ecliptic longitude of Jupiter after some period of 
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time from a selected Earth departure date and that of Earth at the departure date. The 
time T - Tg ^ is the elapsed time after the Earth departure date, and the corresponding 
heliocentric radius to Jupiter is rj. 

The point where the curve for the calculated trajectories intercepts the same depar- 
ture date line in both parts of sketch (a) determines the matching conditions of the Earth 
to Jupiter phase of the desired trajectory. It should be noted that the spacecraft will not 
actually intercept Jupiter, but the orbit of the planet is a convenient boundary point for the 
first part of the heliocentric trajectory calculations. 



Other matching conditions at arrival at Jupiter can be determined by using the 
matched value of r 2 from the right-hand side of sketch (a) in the next four equations, 
which come either from reference 13 or the law of cosines. 


Arrival flight -path angle: 


y 2 = cos 


-1 


r a r p 


l r 2( r a + r p " r 2) 


(Bll) 
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• Arrival heliocentric velocity: 


where Vj^ = — . 

Cl 



Arrival hyperbolic-excess velocity: 


v h,2 2 - V + v 2 2 - 2V J V 2 cos r 2 


(BX2) 

(B13) 


(B14) 


Jupiter to Saturn Trajectories 

Upon arrival at Jupiter, the gravitational field of the planet causes a perturbation 
in the spacecraft trajectory. (See fig. 1.) This perturbation is a strong function of the 
closest passage distance (fig. 2), the value of which must be such that the departing space- 
craft will lie on an intercept trajectory with Saturn. 

The method employed in determining the closest passage distance is somewhat sim- 
ilar to that used for the Earth to Jupiter transfer. A number of trajectories are calcu- 
lated by using the following equations ((B15) to (B29)) with the previously determined 
arrival conditions; the passage distance Rp and the heliocentric radius of Saturn r^ 
are considered as variables. The actual matching procedure used is given after the 
equations. 

The gravity-turn angle with respect to Jupiter, as shown in figure 2 and derived in 
reference 14, is: 


a = 2 sin" 



where V c = 42.581 km/sec for Jupiter. 
For pure gravity turns, 


%2 = V h,3 


(B15) 


(B16) 


because no propulsive velocity is added. 
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The law of cosines can be applied to determine the departure heliocentric velocity , 
and the flight-path angle at Jupiter. The resultant departure heliocentric velocity at 
Jupiter is: 


v 3 2 = V J 2 + v h,3 2 ‘ 2v j v h,3 cos <02 + “) 


(B17) 


= 1 +1 


'h t 3 


'h,3 


- 2 —J— cos (0 2 + a) 
J / V J 


(Bl7a) 


where 


^2 = cos 


■1 PJ 


2 v 2\ 
'h,2 ~ V 2 


2V J V h,2 


V 2 -i-V 2 V 2 ^ 
-lf V J + V 3 " V h,3 

2V J V 3 ' 


(B18) 


and the departure flight -path angle is 

yg = cos 

/V 3 \ 2 

If in equation (B17a), | — | < 2, elliptical equations (B19) to (B22) should be used, and if 

/ x2 V J/ 

/v 3 y 

1 — 1 > 2, hyperbolic equations (B23) to (B29) should be used. In the elliptical and hyper- 
bolic relationships which follow, the matching value of r 2 determined in the right- 
hand side of sketch (a) is used for r^. 

The following elliptical equations are obtained from reference 13: 

Semimajor axis: 


a = 


M 


(B19) 
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* Eccentricity: 


e 2 = 1 — | (2a - rg) cos 2 73 (B20) 

a 2 

For the orbital central angle, eccentric anomaly, and time from periapsis passage, the 
general elliptical equations (B6), (B7), and (B8) are applicable and thus are employed. 

The reference axis of these equations can then be transferred to the departure point at 
Jupiter. The Jupiter-Saturn central angle, obtained by adding or subtracting the central 
angle between the periapsis and the departure point at Jupiter 6 g to or from the central 
angle between the periapsis and the arrival point at Saturn, 8 4 is: 

63 4 = $4 ^3 for TygJ (B21) 

Similarly, the Jupiter to Saturn trip time is 

t 3 4= t 4* t 3 (±t 3 for ?y 3 ) (B22) 

The following hyperbolic equations are obtained from reference 13: 

Semimajor axis: 


a = 


w 


(B23) 


Eccentricity: 


e 2 = 1 + ^|(2a + r 3 ) cos 2 y 3 


(B24) 
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Orbital central angle: 


6 = cos' 


- 1 ( -a + ae 


2 _ 


er 


Hyperbolic anomaly: 

F = cosh-l(-L±£) 


Time from periapsis: 


T = 



+ e sinh 


F) 


(B25) 


(B26) 


(B27) 


As was done for the elliptical transfers, the reference point of the central angle 
and the time of the hyperbolic trajectories between Jupiter and Saturn is transferred to 
the departure point at Jupiter from the periapsis. Thus the Jupiter-Saturn central angle 
is: 


6 3, 4 ~ e 4 ± 6 3 ( ±0 3 for Ty 3) ( B28 ^ 

and the Jupiter-Saturn trip time becomes 

t 3,4 =t 4 ±t 3 (^3 for Ty 3) < B29 ) 

In order to determine the matching trajectory-planetary conditions between Jupiter 
and Saturn, the results from equations (B21) and (B22) or (B28) and (B29) and the ephem- 
eris data (refs. 6 and 7) are plotted as shown in sketch (b). In this sketch the resultant 
angle Zg - Zj 3 represents the difference between the heliocentric ecliptic longitude of 
Saturn after the known departure date from Jupiter and that of Jupiter at the time of depar- 
ture. The time T - Tj 3 is the elapsed time after the Jupiter departure date, and r^ 
is the corresponding radius to Saturn. 

The points in the right-hand plot of sketch (b) where the curves for the calculated 
trajectories cross the Jupiter departure date line determine the values of 0g 4 that are 
used in the left-hand plot of sketch (b) to form a line between the trajectories. The point 
where this line and the departure date line intersect in the left-hand plot determines the 
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matching conditions for 0 3 4 and t 3 4 for the trajectory considered. Then the line 
* connecting the trajectories is replotted so that the passage distance can be obtained 

directly for matching conditions. This procedure completes the determination of all the 
matching conditions of the Jupiter -Saturn phase of the desired trajectory. Now all condi- 
tions have been established for a coplanar trajectory going from Earth to close passage 
of Jupiter and continuing on to Saturn. 
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EQUATIONS USED FOR THE NONCOPLANAR TRAJECTORIES 


Location of Plane Changes 

Additional velocity is required for the coplanar trajectories when the orbital incli- 
nation angle of Jupiter is included. The point of application of this velocity may be con- 
sidered to be either at Earth upon departure or at Jupiter's ascending node. In either 
case, all the matched characteristics from equations (Bl) to (B12) for the coplanar tra- 
jectories between Earth and Jupiter can be maintained by rotating the velocity vector in a 
plane normal to the ecliptic. The geometry used to describe these plane-change maneu- 
vers is presented in figure 3. 


Plane change at Earth .- The inclination angle ij g required for a plane change at 
Earth can be related to the Earth-Jupiter central angle ^ anc * the celestial latitude 
£j of Jupiter at arrival: 


i 


1,E 



(Cl) 


If the heliocentric velocity Vj is held constant, the velocity increment that must be 
added to the coplanar hyperbolic -excess velocity at Earth for a plane change is 


where 


AV h ,i 






cos 


h,E “ V h,l 


(C2) 


and 

Yl = 0 

Plane change at Jupiter's node .- Even when it is desirable to wait until the space- 
craft reaches Jupiter’s ascending node to make a plane change, the trajectory character- 
istics such as semimajor axis a, eccentricity e, and launch date are the same as for the 
coplanar trajectory. The central angle between Earth and the ascending node of Jupiter is 


6 l,n = n j 


± l 


E,1 


(C3) 


where is the longitude of Jupiter's ascending node and lg j is the longitude of 
Earth at departure. (See fig. 3(b).) The velocity of the spacecraft at the node is 
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v n = 


_ M 1 


+ e cos 6 


a(l - e^) 


l,n) _ M 


(C4) 


This velocity must be rotated through an angle equal to the orbital inclination of Jupiter 
*J,E = 1*306°. The rotation of the velocity vector requires a velocity increment 

AV n = 2V n Cos sin (~^) (C5) 

where y n is the flight -path angle of the spacecraft at the node. 

Jupiter to Saturn Trajectories 

As stated earlier, up to the time of arrival at Jupiter, equations (Bl) to (B12) are 
applicable for the noncoplanar trajectories. However, equations (B13), (B14), (B17), and 
(B18) are not applicable when the inclination angle from equation (Cl) is included. There- 
fore, instead of assuming a 0° arrival declination angle, it can be included and determined 
from: 



This expression is based on the geometry presented in figures 3 and 4 and the out-of- 
ecliptic characteristics of Jupiter. 

If this declination angle and all the elements of Jupiter's elliptical orbit are taken 
into account, the arrival hyperbolic-excess velocity at Jupiter becomes 

V h,2 2 = V J 2 + V 2 2 “ 2V J V 2 ( sin y 2 sin y J + cos y 2 cos ^J cos l 2,j) (C7) 

where the terms in the parentheses resulted from using the direction cosines of the vec- 
tors Vj and Vg shown in figure 4. Jupiter's orbital velocity Vj and flight -path 
angle y j were computed from the well-known expressions: 
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(C 8 ) 


y 3 


= cos 


-1 


r a,J r p,J 


r j( r aJ 


+ r 


p,J 


r j) 


(C9) 


Upon departure from Jupiter, the inclination angle can be approximated by 


h,j = sin_1 


( r J Sin *3 J 


(CIO) 


where values from the coplanar results are utilized for rj and 0 g 4 . The departure 
heliocentric velocity Vg can be solved from the following expressions, which are similar 
to equation (C7): 

V h3 = Vj + Vg - 2VjVg (sin yg sin + cos yg cos yj cos ig jj (Cll) 


r T V h 3 r T r o 9 , H 

= -jTp + V 3 ‘ 2 V J V 3 ( sin ^3 Sin y 3 + cos r 3 cos y cos h,j) 


(Cl la) 


r y « r y z 

If - 2 <2 in equation (Cl la), continue with equations (B19) to (B22), and if ■-■■■- - > 2, 

fi fj. 

continue with equations (B23) to (B29). The departure flight -path angle yg and the helio- 
centric radius to Saturn r . are used as variables in computing several trajectories from 

’ 9 

r T V 3 

equation (Cll) and the corresponding equations noted for the conditions of — z — — . These 

M 

results are plotted in a manner similar to that used for the coplanar trajectories in 

D TJ 

sketch (b) of appendix B, where now yg replaces — Then the matched values of 

R J 

yg, 0 g 4 , and tg 4 are determined. 

The remaining matched conditions can be determined for the final phase of the non- 
coplanar Earth-Jupiter-Saturn trajectory by employing these data in the following expres- 
sions and in equation (B15). 

The angle between the arrival hyperbolic -excess velocity vector and Jupiter’s plane 
is 

= -;--lf V 2 COS y 2 sin ^ 


^2 


(C12) 


h ,2 
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Likewise, the angle between the departure hyperbolic -excess velocity vector and Jupiter's 
plane can be expressed as 


= sin-ifc ^f y 3 sin i 3?J \ 
3 \ V M } 


(Cl 3) 


The gravity-turn angle with respect to Jupiter is determined from 
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Parameter 


Present data 


Results of 
reference 2 


Launch date 

Sept. 4, 1977 

Sept. 1977 

V h j, km/sec 

10.52 

10.52 

Tj, km 

7.86 X 10 8 

7.74 X 10 8 

t l,2->' r 

1.421 

1.374 

g, km /sec 

12.68 

12.7 

R p/ R J 

4.2 

4.0 

r g , km 

14.14 X 10 8 

14.15 X 10 8 

‘3,4>y r 

1.558 

1.561 

V h>4 , km/sec 

17.77 

17.8 


17.8 







TABLE II.- COMPARISON OF COPLANAR WITH NONCOPLANAR TRAJECTORIES 


Parameter 

Sept. 5, 1977 launch 

Nov. 8, 

1979 launch 

Coplanar 

Noncoplanar 

Coplanar 

Noncoplanar 

(% M , km/sec 

10.363 

10.363 

10.363 

10.363 

(a Vi, deg 

0 

0 

0 

0 

<a) ti,2. y 

1.456 

1.456 

1.526 

1.526 

( a ^i ? 2> de e 

143.53 

143.53 

144.70 

144.70 

AV h j, km/sec 

0 

0.0014 

0 

0.0893 

h,E> de & 

0 

0.94 

0 

2.26 

Vj, km/sec 

12.954 

12.850 

12.757 

12.449 

Yj> deg 

0 

2.61 

0 

0.51 

( a Vg, km /sec 

13.144 

13.144 

12.770 

12.770 

( a )y 2 > deg 

54.72 

54.72 

54.78 

54.78 

i 2> j, deg 

0 

1.93 

0 

1.85 

deg 

0 

1.28 

0 

1.18 

v h,2> km/sec 

11.996 

11.429 

11.744 

11.510 

Vg, km/sec 

24.462 

24.004 

14.859 

14.584 

r 3 , deg 

10.91 

10.65 

49.62 

50.08 

i 3) j, deg 

0 

1.23 

0 

2.90 

^3, deg 

0 

2.54 

0 

2.36 


1.596 

1.628 

2.337 

2.285 

e 3,4’ de S 

53,65 

54.05 

26.42 

25.80 

a, deg 

93.85 

98.00 

11.87 

11.08 

R p/Rj 

4.65 

4.51 

114.0 

128.0 


(^Quantities assumed to have the same values for both coplanar 
and noncoplanar trajectories. 
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Figure 1.- Coplanar geometric relations of heliocentric trajectories and planetary configurations. 


Jupiter 



Figure 2.- Geometry of coplanar gravity-turn trajectory with respect to Jupiter. 
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Figure 9.- Relative positions of planets with time and approximate angle requirements of planets at launch time for Jupiter gravity-turn 
trajectories to Saturn. V h i = 10.363 km/sec (34 000 fps) to 12.192 km/sec (40 000 fps); y x = 0°. 











(b) V h J = 10.973 km/sec (36 000 fps). 

Figure 12.- Effect of hyperbolic-excess velocity on profiles of gravity-turn trajectories and their relationship to the planets. Launch year, 1978. 
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Q Jupiter* s node 

Gravity turn 


Launch date 


Figure 13.- Effect of plane-change location and trajectory mode on additional velocity requirement for noncoplanar transfers to Saturn 

V h , = 10.363 km/sec (34 000 fps)- ri = 0°. 
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Figure 14.- Effect of Jupiter gravity-turn trajectory on trip time to Saturn for launches during four consecutive years, 
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Closest passage distance, Jupiter radii Turning 






Figure 16.- Effect of launch date on change in both direction and magnitude of heliocentric velocity vector as spacecraft encounters Jupiter. 
V h j = 10.363 km/sec (34 000 fps); = 0°. Flagged symbols represent noncoplanar results. 
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Arrival hyperbolic-excess velocity at Saturn, V h km/sec 


Launch 



0 2 4 6 8 10 1 

Velocity gained at Jupiter, V 3~ V 2* km / sec 


Figure 17.- Effect of increase in velocity of spacecraft at Jupiter on arrival hyperbolic-excess velocity at Saturn. 

V h ! = 10.363 km/sec (34 000 fps); y, = 0°. 




<b) 1977 launch. 

Figure 18.- Profiles of Earth-Jupiter- Saturn trajectories for 1976 to 1979 launches. V h i = 10.363 km/sec (34 000 fps); Yi = 0°. 
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Figure 20. Effect of flight-path angle at Earth departure on launch window. V. , = 10.363 km/sec (34 000 fps). 
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Figure 23.- Effect of hyperbolic-excess velocity and flight-path angle on trip time. Launch year, 1978. 
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Figure 24.- Comparison of injected payload for a direct mission with a Jupiter gravity-turn mission to Saturn. I sp - 4177.8 m/sec 

(426 Ibf-sec/lbm); o = 0.105. 


NASA- Langley, 1966 L-5011 
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" The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space . The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 
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